Allosteric inhibition of the epidermal growth factor receptor through disruption of transmembrane interactions

The epidermal growth factor receptor (EGFR) is a receptor tyrosine kinase (RTK) commonly targeted for inhibition by anticancer therapeutics. Current therapeutics target EGFR’s kinase domain or extracellular region. However, these types of inhibitors are not specific for tumors over healthy tissue and therefore cause undesirable side effects. Our lab has recently developed a new strategy to regulate RTK activity by designing a peptide that specifically binds to the transmembrane (TM) region of the RTK to allosterically modify kinase activity. These peptides are acidity-responsive, allowing them to preferentially target acidic environments like tumors. We have applied this strategy to EGFR and created the PET1 peptide. We observed that PET1 behaves as a pH-responsive peptide that modulates the configuration of the EGFR TM through a direct interaction. Our data indicated that PET1 inhibits EGFR-mediated cell migration. Finally, we investigated the mechanism of inhibition through molecular dynamics simulations, which showed that PET1 sits between the two EGFR TM helices; this molecular mechanism was additionally supported by AlphaFold-Multimer predictions. We propose that the PET1-induced disruption of native TM interactions disturbs the conformation of the kinase domain in such a way that it inhibits EGFR’s ability to send migratory cell signals. This study is a proof-of-concept that acidity-responsive membrane peptide ligands can be generally applied to RTKs. In addition, PET1 constitutes a viable approach to therapeutically target the TM of EGFR.

The epidermal growth factor receptor (EGFR) is a receptor tyrosine kinase (RTK) commonly targeted for inhibition by anticancer therapeutics. Current therapeutics target EGFR's kinase domain or extracellular region. However, these types of inhibitors are not specific for tumors over healthy tissue and therefore cause undesirable side effects. Our lab has recently developed a new strategy to regulate RTK activity by designing a peptide that specifically binds to the transmembrane (TM) region of the RTK to allosterically modify kinase activity. These peptides are acidity-responsive, allowing them to preferentially target acidic environments like tumors. We have applied this strategy to EGFR and created the PET1 peptide. We observed that PET1 behaves as a pH-responsive peptide that modulates the configuration of the EGFR TM through a direct interaction. Our data indicated that PET1 inhibits EGFR-mediated cell migration. Finally, we investigated the mechanism of inhibition through molecular dynamics simulations, which showed that PET1 sits between the two EGFR TM helices; this molecular mechanism was additionally supported by AlphaFold-Multimer predictions. We propose that the PET1-induced disruption of native TM interactions disturbs the conformation of the kinase domain in such a way that it inhibits EGFR's ability to send migratory cell signals. This study is a proof-ofconcept that acidity-responsive membrane peptide ligands can be generally applied to RTKs. In addition, PET1 constitutes a viable approach to therapeutically target the TM of EGFR.
The epidermal growth factor receptor (EGFR) is a HERfamily receptor tyrosine kinase (RTK) that is involved in cell signaling in healthy tissue. Activation of EGFR regulates essential cellular processes including cell migration, proliferation, and apoptosis (1). To mediate these processes, the extracellular ligand-binding region of EGFR senses environmental cues via interactions with one of its seven known ligands, of which epidermal growth factor (EGF) is the most well characterized (2,3). Ligand binding promotes EGFR oligomerization mediated by the extracellular region. Signaling is then transduced across the membrane by altering the configuration of the transmembrane (TM) domain, by dimerization of the TM helical region or a change in the arrangement of the TM helices within such a dimer. Specifically, the TM of unliganded (inactive) EGFR dimerizes at the C-terminus (C t ), while the ligand-bound form dimerizes N-terminally (N t ), and the two helices are also rotated by 180 between the conformations (4,5). The ligand-bound TM configuration promotes asymmetric dimerization of the intracellular juxta-membrane (JM) and kinase domains, which causes autophosphorylation of intracellular tyrosine residues (6,7). Effector proteins are then recruited to activate various cellular signaling pathways, including RAS/RAF/MEK, PI3K/ AKT/mTOR, and JAK/STAT (1).
Because of its essential roles in cell signaling, misregulation or overexpression of EGFR often causes a cancerous phenotype. Indeed, EGFR is commonly overexpressed in solid tumors, such as breast, colon, head-and-neck, renal, ovarian, and non-small-cell lung cancer (8)(9)(10)(11). Furthermore, EGFRmediated cancers tend to be more aggressive (12). Currently, there are two main therapeutic approaches that are effective for targeting EGFR in cancer: monoclonal antibodies and small-molecule tyrosine kinase inhibitors (11). Both approaches are generally safer and more efficacious than chemotherapy. However, major challenges for antibodies include a short lifespan and variations in tumor development (13), and tyrosine kinase inhibitors are often promiscuous among other RTKs due to the highly conserved ATP-binding pocket causing off-target effects (11,14). Additionally, both strategies often become less effective over time as tumors develop resistance (11,13,14). For these reasons, it is necessary to find safer, more effective, and more selective ways to inhibit EGFR activity.
Recently, our group has developed a novel approach to modify RTK activity by targeting the receptor's TM domain using a pH-responsive peptide (15). Acidity-responsive peptides such as pHLIP (16) and ATRAM (17,18) are marginally hydrophobic and contain acidic residues across the sequence.
At physiological pH, the acidic residues are unprotonated and therefore negatively charged, allowing the peptide to be soluble, but able to bind to the surface of lipid membranes, in an unstructured conformation. At lower pH, the acidic residues become protonated, resulting in membrane insertion and a gain of α-helical structure. These peptides can preferentially target cancer cells over healthy tissue by taking advantage of the slightly acidic extracellular pH that is a hallmark of tumors (18)(19)(20)(21)(22). Alves et al. (15) evolved this concept and designed the peptide TYPE7 to be specific for the RTK EphA2. TYPE7 binds the TM to allosterically regulate EphA2 kinase activity by causing a configurational change (23). This method of regulation is likely to be more selective for EphA2 than targeting the highly conserved kinase domain. When combined with the increased selectivity for cancer cells, TYPE7 represents a potentially useful development to target RTK activity in cancer.
We sought to use the TYPE7 approach to inhibit EGFR. Here we report a novel pH-responsive Peptide for EGFR Targeting (PET1). PET1 binds selectively to EGFR in cancer cells and inhibits the ability of the ligand EGF to activate cell migration. Interestingly, PET1 does not modify the receptor's oligomerization state. Molecular dynamics simulations and AlphaFold-Multimer predictions reveal that PET1 disrupts EGFR by forcing apart the TM dimer, bridging the individual TM domains. PET1 therefore induces a configuration of EGFR that is inactive without full dissociation of the oligomeric complex.

Results
PET1 is a pH-responsive peptide that interacts with the TM region of EGFR To design PET1, we modified the TM region of human EGFR as previously described (Fig. 1A) (15). Glutamic acid (E) residues were strategically placed throughout the TM region and at the charged JM region immediately C t to the TM. Glutamic acid residues are negatively charged at neutral pH and only become protonated at low pH and therefore are expected to confer pH-responsiveness to the peptide. To Figure 1. PET1 is a pH-responsive transmembrane peptide that binds to the TM of EGFR. A, the sequence of EGFR (residue 638-683) with the TM domain underlined is aligned with the sequence of PET1 and TM-EGFR peptides. Amino acid mutations to acidic residues are highlighted in red. The C t CWN tag on TM-EGFR is highlighted in blue. B, CD spectra informs on the secondary structure changes of PET1 in buffer at pH 7.5 (black), in POPC vesicles at pH 7.5 (gray), and in POPC vesicles at pH 4.2 (red). C, oriented circular dichroism of PET1 in POPC supported bilayers at pH 4.2 (red). D, the fluorescence center of mass of PET1-NBD in POPC vesicles alone (black, top cartoon) or containing TM-EGFR (pink, bottom cartoon) was determined at varying pH values to determine the pH 50 of insertion (dashed lines) using Equation 2. Reported pH 50 values (mean ± S.D.) are an average of three individual replicates. Statistical analysis was performed using a t test (p = 0.029). E, the fluorescence spectra was recorded for the C t W residue of TM-EGFR in POPC vesicles at pH 4.3 and 7.4 in the presence (pink) or absence (gray) of PET1. Box plot conveys the wavelength corresponding to the maximum fluorescence of the curve. N = 6. Statistical analysis was performed using a one-Way ANOVA (p = 3 × 10 −5 ) with a post hoc Dunnet-T3 for comparisons between groups, as the Levene Statistic was significant (p = 0.015). EGFR, epidermal growth factor receptor; PET1, peptide for EGFR targeting; POPC, 1-palmitoyl-2-oleoyl-sn-glycero-3phosphocholine; TM, transmembrane. determine the interaction of PET1 with lipid membranes, we performed complementary biophysical assays. We used CD and oriented CD (OCD) to assess PET1 secondary structure and the average tilt of membrane insertion, respectively (24). CD performed in buffer (Fig. 1B, black) and 1-palmitoyl-2oleoyl-sn-glycero-3-phosphocholine (POPC) lipid vesicles at pH 7.5 (Fig. 1B, gray) revealed PET1 is in a random coil conformation in both conditions, as indicated by the minima at 200 nm. However, in vesicles at pH 4.2, PET1 adopts an αhelical conformation based on the minima at 208 and 222 nm (Fig. 1B, red), indicating that PET1 is only structured at acidic pH. To determine if folding is due to membrane insertion, we performed OCD at acidic pH. OCD measures the average tilt relative to the membrane normal of an α-helix. The OCD curve of PET1 (Fig. 1C) is consistent with a peptide configuration that is inserted and oriented with a noticeable helical tilt, as indicated by the similar intensity at 208 and 225 nm. Such orientation is consistent with NMR structures of the EGFR TM domain (4,25,26). However, the OCD method does not allow to distinguish between a stable peptide tilt and the presence of two conformations of different membrane orientation, which might be dynamically interconverting. Together, the biophysical results reveal that PET1 is pH-responsive because the peptide undergoes the desired shift from unstructured to TM when the pH decreases.
We performed a pH titration experiment to further investigate the pH-dependent membrane insertion of PET1. A fluorescently labeled PET1 (PET1-NBD, Fig. S1) was incubated with POPC vesicles, and the pH was changed to cover a wide pH range from acidic to neutral values. NBD is an environmentally sensitive dye that presents a blue-shifted fluorescence spectra in a hydrophobic environment (27)(28)(29). Therefore, we used the NBD fluorescence center of mass (COM, Equation 1) as an indicator for the peptide environment ( Fig. 1D). At low pH, we observed a low COM that suggests a more hydrophobic (probably membraneassociated) environment, while the higher COM at neutral pH suggests that the NBD in PET1 is more solvent-exposed. The transition between states occurred in a sigmoidal fashion with a midpoint (pH 50 ) of 5.00 ± 0.01 (24). For pH-responsive peptides that bind the TM of an RTK, the presence of that RTK can increase the pH 50 due to an increase in tendency to be inserted in the presence of a membrane-binding partner (15,23). For this reason, we repeated the experiment using proteo-liposomes containing a peptide mimic of the TM of EGFR (TM-EGFR). Under these conditions, we observed that the pH 50 value increased to 5.36 ± 0.16 (Fig. 1D). This result indicates that the presence of TM-EGFR makes PET1 insertion more favorable, suggesting that PET1 interacts with the TM region of EGFR. PET1-NBD experiments described the effect that TM-EGFR causes in PET1 insertion. We performed an orthogonal experiment to determine if PET1 also affects the configuration of TM-EGFR. For this, we used the C t tryptophan (W) residue on TM-EGFR as a fluorescent reporter of hydrophobicity, similarly to NBD. W also presents a blue shift in more hydrophobic environments (30). Using pH values representative of the fully TM (pH 4.3) or fully unstructured (pH 7.4) PET1 baselines as determined by the titration experiment, we measured the W fluorescence spectra of TM-EGFR in POPC lipid vesicles with and without PET1 added. We used the spectral maximum as an indicator for W positioning. At pH 4.3 when PET1 is fully inserted, the addition of PET1 caused a significant increase in the spectral max wavelength (Fig. 1E). This wavelength increase was accompanied by a significant fluorescence decrease (Fig. S2). The observed spectral red-shift indicates a transition to a more polar environment when PET1 forms a TM helix. At pH 7.4 when PET1 does not insert to the membrane, the addition of PET1 had no effect (Figs. 1D and S2). These results suggest that only the TM conformation of PET1 modifies the environment of the TM-EGFR C t , due to an interaction between PET1 and TM-EGFR.

PET1 inhibits EGFR-mediated cell migration
We next sought to determine if PET1 modifies EGFR activity. We used cell migration as an indicator of downstream EGFR-regulated cell signaling, since activation of this RTK promotes cell migration. We performed a Boyden cell chamber assay in which A375 melanoma cells migrate through a porous membrane in response to a chemoattractant (Fig. 2). As expected, treatment with EGF significantly enhanced cell migration due to its ability to activate EGFR (31,32). We observed that PET1 alone did not change basal levels of migration, but interestingly, PET1 was able to significantly block the ability of EGF to promote migration. We have previously shown that pHLIP, a peptide with similar pH-responsive properties to PET1, is not able to affect cell migration of A375 cells (15), which suggests that effect of PET1 is specific to its interaction with EGFR. To validate the migration results, we performed a control cell viability assay using the 3-(4,5-Dimethylthiazol-2yl)-5-(3-carboxymethoxyphenyl)-2-(4-sulfophenyl)-2H-tetrazolium (MTS) reagent (Fig. S3). We found that PET1 caused no significant cellular toxicity, with or without EGF treatment. Our data indicate that PET1 inhibits ligand-induced activation of EGFR.

PET1 does not alter EGFR oligomerization
We investigated next the mechanism through which PET1 inhibits EGFR activation. Activation of EGFR by EGF promotes receptor self-assembly by stabilizing the dimeric and oligomeric states. We therefore studied the effect of PET1 on EGFR oligomerization via pulsed interleaved excitation fluorescence cross-correlation spectroscopy (PIE-FCCS) (33). PIE-FCCS is a time-resolved fluorescence method in which two excitation lasers are focused on the plasma membrane of live cells, and fluorescence fluctuations are recorded to quantify the expression level, mobility, and oligomerization state of the labeled membrane proteins. Single cell data are fit to determine the fraction of the codiffusing species relative to the red or green species (f c ). We assessed the oligomerization state of EGFR before and after PET1 addition by comparing f c values (Fig. 3A). In order to determine the effect of the peptide on EGFR receptor-receptor interactions, we performed two sets of experiments. In the first, we collected data after PET1 was added to EGFR-expressing cells in the absence of ligand. Then, we added EGF to determine if the peptide affected ligandstimulated multimerization of EGFR. In unstimulated cells, the median f c value was 0.00, indicating that EGFR was predominantly monomeric. Upon peptide addition, the median f c value was 0.03, with no statistical difference compared to unstimulated EGFR. Upon EGF addition to these PET1-treated cells, a median f c value of 0.21 was obtained. This value is significantly larger, confirming that EGF considerably promotes ligandstimulated EGFR multimerization even in the presence of PET1.
The second set of the experiments was to test if a multimeric state, generated by first adding EGF ligand to EGFR, could be disrupted by PET1. First, we collected data after EGF was added to EGFR-expressing cells. We then added PET1 to the media and collected more PIE-FCCS measurements. EGF-stimulated EGFR yielded a median f c value of 0.19. When the peptide was added to the well, the median f c values were 0.18. There was no significant difference between the f c values (Fig. 3A) of EGFR oligomers, regardless of the presence of the peptide, indicating that PET1 does not disrupt the oligomerization of EGFR. In fact the diffusion of the EGFR:PET1 complex is slightly slower than the diffusion of the EGFR oligomer alone (Fig. 3B).

PET1 colocalizes with and binds to EGFR
To demonstrate that the effect of PET1 was direct, we sought to demonstrate the interaction between PET1 and EGFR. We first determined the cellular location of PET1 with respect to EGFR by treating A431 cells, which have a naturally high level of EGFR expression, with a fluorescently tagged  version of PET1 (PET1-DL680, Fig. S1). We visualized the cellular location of PET1 and EGFR via confocal microscopy ( Fig. 4A). PET1 localized at the plasma membrane, and there was a strong overlap between PET1 and EGFR signals irrespective of EGF treatment. To quantify colocalization, we determined the Pearson's correlation coefficient (r) (Fig. 4B). The r parameter ranges from +1 for perfect correlation to −1 for anticorrelation (34). The value of r was 0.7 with and without EGF, indicating a strong correlation between PET1-DL680 and EGFR cellular localization regardless of whether EGF is present. This result suggests that PET1 can bind to EGFR prior to engagement with EGF.
Next, we performed co-immunoprecipitation (co-IP) experiments to investigate whether colocalization was due to EGFR-PET1 binding. We immunoprecipitated EGFR from cells treated with PET1-DL680 using the gentle detergent NP-40 to allow for precipitation of bound proteins. Figure 4C shows that PET1-DL680 at varying concentrations co-IP's with EGFR. We also observed that binding is independent of the presence of EGF, in agreement with the colocalization results. In addition to the low molecular weight band corresponding to PET1-DL680 (broad band at 4 kDa), we observed a high molecular weight fluorescent band likely due to an SDSresistant complex of PET1-DL680 and EGFR (180 kDa). Co-IP and SDS-resistant binding suggest that the interaction between PET1-DL680 and EGFR is strong.
To determine if PET1 action on EGFR is selective, we utilized an array of 58 RTKs to measure phosphorylation in control conditions and in the presence of EGF. We studied the effect of PET1 and used a scrambled peptide (SP) as a negative control (Fig. S4). PET1 alone did not cause phosphorylation of any RTKs, suggesting that the peptide does not cause off-target effects. We observed that after EGF treatment, EGFR was heavily phosphorylated (red) as expected, and EGFR coreceptors showed increased phosphorylation including ErbB2/ Her2 (orange), ErbB3/HER3 (purple), MerTK (green), and EphB2 (yellow) (35)(36)(37). Interestingly, while EGFR phosphorylation was not sensitive to the presence of the peptide, phosphorylation of MerTK and EphB2 appeared to show a decrease with EGFR+PET1 treatment. EGFR-mediated phosphorylation of coreceptors may be one reason for changes in downstream signaling in the latter cases.

PET1 separates the EGFR TM dimer by binding both helices simultaneously
Once we determined that PET1's effect on cell migration was due to binding between PET1 and EGFR, we sought to further understand how PET1 binds the EGFR TM through molecular dynamics simulations (Fig. 5). We employed Coarse Grain (CG) molecular simulations for our study to observe the binding of PET1 with the EGFR TMs. For this, we studied the . PET1 interacts with endogenous EGFR. A, representative images of A431 cells treated with PET1-DL680 (magenta) for 1 h followed by a 5 min incubation with or without EGF, fixed, and stained for EGFR (green). DAPI was used for nucleus staining (blue). Scale bar represents 75 μm. B, Pearson's correlation coefficient (r) was calculated for PET1 and EGFR channels. N = 3, n = 15. The error bars represent SEM. C, EGFR was immunoprecipitated from lysates of A431 cells treated with PET1-DL680 for 1 h followed by 5 min treatment with or without EGF. Control (CT) lane shows that PET1-DL680 alone runs as a wide band. SDS-PAGE and Western blot of the eluates for EGFR (green) is shown. EGF, epidermal growth factor; EGFR, epidermal growth factor receptor; PET1, peptide for EGFR targeting.
Inhibition of EGFR via a transmembrane peptide association of EGFR-TM regions in the absence and presence of PET1 (Fig. 5, A and B). As the starting configurations, we modeled the membrane-embedded regions as ideal monomeric helices placed in random orientation 5 nm apart from each other in a POPC membrane bilayer. The CG simulation was run for 4 μs in quadruplicates to assess the consistency of the results. In the case of the EGFR TM alone (Fig. 5C), we saw that the two helices came together within 0.5 μs. Table S1 shows the PREDDIMER analysis for the EGFR TM dimers in comparison with the NMR structure (PDB ID: 5LV6) (5). In presence of PET1, the Fscor-a score reporting on the close and energetically strong packing of helices-was very low for the most populated cluster centers (two with 0, the two lesser populated structures with 3.5 and 2.6, respectively; Table S1), suggesting weak association between the EGFR TMs for most of the simulations. However, in the control simulations without PET1, the Fscor was high (3.2-4.0), suggesting stable interaction between the EGFR TMs. In the absence of PET1, the EGFR TMs assume a configuration which on average resembles that of the NMR structure thought to be the active state (RMSD of top cluster 1.7 Å, Table S1). This result, as well as CG simulations for other single membrane crossing TM receptors, gives us confidence that the employed Martini 3.0 potential function leads to accurate EGFR TM dimer structures.
For the simulations with PET1, we assumed the simplest possible scenario, where PET1 binds to the EGFR TM dimer. We observed the association of the PET1 monomer with the EGFR TM helices forming a heterotrimer within 1 μs in the four molecular dynamics trajectories (Fig. 5D). To better observe the effect of PET1 association, we RMSD-clustered the simulations and superimposed the main conformer in Figure 6, A and B. We also show the contact maps in Figure 6, C and D, as averages over all simulations and plotted as the interhelical distance between the EGFR TMs (considering the COM of the TM region) versus the crossing angle between the two EGFR TM helices in Figure 6, E and F, yielding a 2D population map. In all cases, the helices aligned close to parallel. In the EGFR TM dimer structure, the COM of the helices was spaced 0.65 to 0.95 nm, while with PET1, the distance increased to 0.90 to 1.15 Å or 1.30 to 1.55 nm. Thus, overall the presence of PET1 prevents the direct association of the EGFR TM helices for most of the time in the simulations as shown in Figure 6B. This is true in three of the four cluster centers, but in the case of the third cluster, center PET1 has a tilted orientation in the bilayer and associates more on the side of a weak EGFR TM dimer (Fig. 6B), leading to the higher Fscor, as mentioned above. We also plotted the configurational transition of EGFR TM dimers in the absence and presence of PET1 (Fig. S5), demonstrating that both systems are somewhat dynamic, especially the EGFR TM configurations in the presence of PET1 which undergoes considerable fluctuations, most of which are accompanied by temporary increases in helix-to-helix COM distances.
As a control, we repeated the CG simulations with the SP PET1 (Fig. S4). The CG results showed that the SP was able to bind to the EGFR TMs (Fig. S6), but this interaction did not robustly separate the EGFR dimer. In the presence of SP, the TM helices largely remained within 0.75 nm (Fig. S6C), similarly to the case when the simulations were run in the absence of peptide (Fig. 6E). Additionally, SP did not prevent key residues in the EGFR TMs to engage in dimer-stabilizing interactions (Fig. S6B). These simulations suggest that the effect of PET1 is specific. To benchmark the CG simulations, we used AlphaFold-Multimer (38) to predict the structure of the complex formed by the EGFR TM helices and PET1 in a 2:1 stoichiometry. The artificial intelligence program predicted that PET1 binds to both TMs and blocks their C t association (Fig. S7), in agreement with the CG simulations. As a control, we also applied AlphaFold-Multimer to the isolated EGFR TM helices. The obtained prediction is in strong agreement with the TM structure believed to correspond to the active conformation of the receptor, solved by NMR (Fig. S7) (26). The confidence score of both predictions had reasonable values (0.58 versus 0.43). However, the AlphaFold-Multimer prediction with SP yielded a low score (0.3) and is therefore not considered a robust prediction. To summarize, AlphaFold-Multimer and the CG simulations support the notion that PET1 disrupts native interactions between the TM helices of EGFR.

Discussion
Our data indicate that PET1 is a pH-responsive peptide that inhibits EGFR activation through disruption of the EGFR TM dimer. The pH-responsive nature of PET1 is an important aspect of its design. Acidity-responsive peptides like pHLIP can preferentially target tumors over healthy tissue in mouse models, due to the more acidic environment of cancer cells (18)(19)(20)(21)(22). While EGFR is commonly overexpressed in cancer, there is still expression in healthy tissue. Thus, designing PET1 as a pH-responsive peptide will likely ameliorate any eventual off-target effects. We observed that PET1 displayed an increased pH 50 in the presence of TM-EGFR (Fig. 1D). However, the reported pH 50 of PET1 in the presence of TM-EGFR is lower than the extracellular pH of tumors (pH 6.4-6.8). One must consider that pH 50 determination by an N t NBD tag consistently provides a lower pH 50 value than determination by tryptophan spectral max, COM, or CD (24). Using this difference of 0.5, we can correct the determined pH 50 in the presence of TM-EGFR to what is likely a more accurate pH 50 of 5.85, which is closer to that of cancer cells. In addition, the pH 50 was determined in pure POPC lipid vesicles. It is known that the lipid composition of a membrane can strongly affect the pH 50 of pH-responsive peptides (18). Therefore, this Inhibition of EGFR via a transmembrane peptide in vitro model membrane may not directly recapitulate the insertion conditions of the peptide into mammalian cell membranes (39). Regardless, our cellular results show that PET1 is able to insert into the plasma membrane and interact with EGFR at physiological pH, similarly to results obtained for the TYPE7 for the modulation of the EphA2 receptor (15).
Our model of PET1 inhibition of EGFR is shown in Figure 7. In its ligand-free form, EGFR primarily exists in an equilibrium between a monomer and dimer. Ligand binding shifts the equilibrium towards the dimer and induces a configurational change that allows kinase activation. In the inactive dimeric configuration, NMR reveals that the EGFR TM dimer has a helix-helix crossing angle of 30 utilizing the C t AxxxG motif (PDB ID: 2M0B) (5). This configuration causes the C t ends of the TM to be positioned only 7.2 Å apart (4,40,41). The JM is therefore held too close together to form the antiparallel dimer characteristic of the active form, and the positively charged residues of the JM interact with the negatively charged lipids of the inner membrane leaflet (26,39,41). These electrostatic interactions hold the kinase domain close to the membrane where it is inactive (26). In contrast, binding of a ligand such as EGF causes the TM dimer to shift positions such that it dimerizes through the SxxxGxxxA motif at the N t of the TM domain, with a helical crossing angle of −42 (4, 5, 26). Moreover, the N t motif is on the opposite side of the helices, compared to the C t AxxxG motif, which is utilized in the inactive state, thus leading to a rotation of both helices by 180 for receptor activation (Fig. 7). In the active dimer, the C t end of the helices are 20 Å apart, a sufficient distance to allow JM antiparallel dimerization, which promotes asymmetric dimerization and activation of the kinase domain (39,40). Our MD data in the presence of PET1 reveals a TM configuration unlike either the active or inactive state; PET1 is sandwiched between the two EGFR helices, which disallows most intramolecular contacts between the EGFR TM regions, and forces the COM of the TM regions to be approximately 15 Å apart. Furthermore, in this configuration, the negatively charged amino acids just outside the membrane on the PET1 sequence likely interact with the positive residues in the JM region of EGFR (not examined here by modeling, since only the very N t region of the JM was included in the simulations). This is supported by our tryptophan fluorescence data (Fig. 1E), which shows the EGFR JM residues in a more solvent-exposed position in the presence of PET1. It is possible that PET1 forces the JM away from the membrane, due to an electrostatic attraction between the acidic residues of PET1 and the basic residues of the JM. This interaction might force a conformation that is neither close enough to force the JMs apart and to bind with the membrane nor far enough to allow the JM to dimerize. We therefore propose that PET1 promotes a conformation that is different to the ligand bound or unbound dimer and therefore disallows the native downstream effects of EGFR. Indeed, the model proposed here is supported by previous findings. Prior MD simulations revealed that the EGFR TM is likely able to exist in configurations other than the active or inactive ones discussed above (42,43). Additionally, cryo-EM has shown that another EGFR ligand, TGF-α, induces an extracellular conformation that is different to the EGF-bound state and likely causes an intermediate TM conformation somewhere between the two previously discussed (44). This is further supported by crosslinking experiments in live cells that reveal that the EGFR TM-JM region is configured differently by binding of each of the seven known ligands (43). Together, these experiments suggest that the TM dimer is more dynamic than originally thought, and it is further proposed that these Figure 7. Model of the EGFR configurational changes caused by EGF and PET1. It is well established that EGFR exists in a monomer dimer equilibrium in its inactive state. The addition of EGF causes a configurational rearrangement for the extracellular, TM, JM, and kinase domains that allow autophosphorylation and activation. We propose that PET1 induces a configuration of the protein in which the TM interaction modifies the ability of the JM and kinase domain to arrange correctly for signaling to occur. EGF, epidermal growth factor; EGFR, epidermal growth factor receptor; JM, juxta-membrane; PET1, peptide for EGFR targeting; TM, transmembrane. small dynamic changes in the TM dimer exert large effects on the intracellular configurations and signaling of EGFR (43). Therefore, it is not at all unreasonable to expect that disruption of TM dimerization has a strong effect on EGFRs activity.
Intriguingly, PET1's mechanism of activity is opposite to that of TYPE7, the only other pH-responsive peptide published to date that targets an RTK (EphA2) (15). TYPE7 works by "stapling together" both TM helices of the EphA2 dimer, which stabilizes the ligand-bound conformation and promotes downstream EphA2 signals (15). In contrast, PET1 disrupts TM binding to inhibit downstream EGFR signals. The opposing mechanism is interesting, as the peptides were similarly designed, with the E mutations placed on the helix interface that participates in ligand-independent dimerization. In contrast, the EphA2 TM homologous peptide N3, a variant of TYPE7 with the E residues placed on the interface that participates in ligand-dependent signaling, appears to function more similarly to PET1 (23). N3 also sits between the two helices and disrupts the EphA2 dimer entirely. PET1, TYPE7, and E3 each contain a tail of acidic residues that can form stabilizing interactions with the positively charged JM region of the targeted receptor, but it seems it is the placement of residues within the TM helix which confers specificity of mechanism. As discussed above, it appears that disruption of the dimer only affects the ability of the RTK to be ligandactivated, not the basal levels of activity. Further work will be needed to understand how to fine-tune the design, so as the peptide stabilizes or disrupts a specific dimer conformation.
Our work is a proof-of-principle study that shows that targeting the TM of EGFR can lead to an efficient inhibition of this receptor. We describe the molecular mechanisms through which PET1 functions, in which PET1 biases the dimerization of the EGFR TM domain to allosterically regulate downstream function. Future work would benefit from fully characterizing the cellular mechanism of PET1: what phosphorylation patterns and signaling pathways are affected? Are other cell phenotypes besides migration affected? In addition, optimization of the peptide's ideal concentration, half-life, and kinetics would be invaluable for better understanding PET1's function. Finally, it would be particularly interesting to determine the stoichiometry of the peptide and EGFR TM complex, which is currently unknown. With further exploration of PET1, we might find that we have a new way to therapeutically target EGFR-mediated cancers that combats off-target effects and drug resistance.

Peptide dye conjugation
For NBD and DL680 conjugation of PET1, the esterified version of each dye (Thermo Fisher Scientific) was linked to the N-terminus. Dye suspended in dimethylformamide (DMF) was added to PET1 dissolved in 100 mM sodium phosphate, 150 mM sodium chloride (pH 7.0) at a dye to peptide molar ratio of approximately 1:10 for NBD and 1:5 for DL680. The mixture was shaken for 1.5 h and then centrifuged at 14,000g to remove precipitated dye. The supernatant was then run on a PD10 desalting column with 1 mM sodium phosphate (pH 7.2) buffer to separate free dye from conjugated peptide.

Liposome preparation
Lipids were purchased from Avanti Polar Lipids. POPC stocks were suspended in chloroform. Aliquots were dried under a stream of argon gas and then subjected to vacuum at least 2 h before resuspension in 10 mM sodium phosphate (pH 7.4). For proteo-liposomes containing TM-EGFR, stocks of TM-EGFR in methanol were mixed with POPC prior to drying. Drying was performed in 13 mm glass culture tubes that had been piranha (75% H 2 SO 4 , 25% H 2 O 2 ) cleaned for 3 min to reduce peptide sticking to the glass. Resuspended samples were extruded using a Mini-Extruder (Avanti Polar Lipids) through a 100 nm membrane (Whatman) to form large unilamellar vesicles (LUVs).

Circular dichroism
Stocks of POPC and PET1 were prepared in chloroform and 1 mM sodium phosphate (pH 7.4) buffer, respectively. An aliquot of POPC was dried under a stream of argon gas before placed in a desiccator for at least 2 h. The POPC film was resuspended with 1 ml of 1 mM NaPi pH 7.4 buffer and extruded through a 100 nm Nuclepore Track-Etch Membrane (Whatman) to produce LUVs. PET1 was diluted to a working concentration of 7 μM peptide suspended in 20 mM sodium phosphate pH 7.5 or 20 mM sodium acetate pH 4.3. PET1 was incubated with LUVs at a 150:1 lipid to Inhibition of EGFR via a transmembrane peptide peptide molar ratio. Samples were recorded on a Jasco J-815 CD spectrometer using a 2 mm quartz cuvette (Starna Cells Inc). All conditions were averaged over two technical replicates. Appropriate buffer backgrounds were collected on the same day and subtracted appropriately.

Oriented circular dichroism
Stocks of POPC and PET1 were suspended in chloroform and 1,1,1,3,3,3-Hexafluoro-2-propanol (HFIP), respectively. An aliquot of POPC was dried under a stream of argon gas before placed in a desiccator for at least 2 h. The POPC film was resuspended with a calculated volume of PET1 stock solution to reach a 50:1 lipid to peptide molar ratio and dried correspondingly. The POPC-PET1 film was resuspended with HFIP and deposited homogenously across two circular quartz slides (Hellma Analytics) cleaned with piranha solution. These slides were placed in glass petri dishes and balanced horizontally within a chemical hood at room temperature overnight to ensure complete HFIP evaporation. Lipid films on each slide were rehydrated with 150 ml of 100 mM sodium acetate buffer pH 4.24 for 16 h in a 96% relative humidity chamber packed with saturated K2SO4. The majority of buffer was removed, and the slides were assembled into an OCD cell packed with saturated K 2 SO 4 to maintain humidity. The OCD spectra were recorded on a Jasco J-815 CD spectrometer and averaged over eight 45 rotations of the cell. POPC lipid backgrounds were collected separately and subtracted appropriately.

pH titration assay
LUVs and proteoliposomes prepared as above in 10 mM sodium phosphate (pH 7.4) were incubated for at least an hour with PET1-NBD at a lipid:TM-EGFR:PET1 molar ratio of 1000:5:1. Stocks were then diluted into a series of 100 mM sodium phosphate or sodium acetate buffers at pH's between 4 and 7.6 in 0.2 intervals. The NBD fluorescence spectra were recorded at 25 C with excitation at 470 nm and an emission range of 520 to 600 nm using a Cytation 5 imaging plate reader (Biotek Instruments). Lipid blanks were prepared at the highest and lowest pH, averaged, and subtracted from the test spectra. The fluorescence (I i ) and wavelength (λ) of each curve was used to calculate the COM at each pH using Equation 1.
The COM was plotted against pH to determine the pH 50 using Equation 2.
1þ10 mðpH−pH 50 Þ Tryptophan fluorescence assay LUVs and proteoliposomes were prepared as above. Appropriately, pH adjusted 100 mM sodium phosphate (pH 7.4) or sodium acetate (pH 4.3) buffer and PET1 were added to LUVs for a final concentration of 200 μM POPC, 1 μM TM-EGFR, and 5 μM PET1. Samples were incubated for 1 h at room temperature (19-21 C) to allow peptide binding to come to equilibrium. Tryptophan fluorescence spectra were then obtained on a Cary Eclipse Fluorescence Spectrophotometer at an excitation wavelength of 280 nm (Agilent Scientific). For all treatments lipid blanks were subtracted.
Cell culture A375, A431, and Cos7 cell lines were obtained from ATCC and maintained at 5% CO 2 and 37 C in Dulbecco's Modified Eagle's Medium (DMEM) supplemented with glucose, 10% fetal bovine serum (FBS), and 100 U/ml penicillinstreptomycin. Cells were passed at 80% confluency and were not used beyond 40 passes. Cell lines were STR tested for authentication via ATCC.
Colocalization A431 cells were plated at 80% confluency on a #1.5 glass coverslip, allowed to adhere for 24 h, and then starved overnight. Cells were then treated with serum-free DMEM without (no treatment) or with (PET1) PET1-DL680 for 1 h prior to a 5 min EGF 100 ng/ml treatment. Cells were washed with PBS containing 1 mM MgCl 2 and 100 mM CaCl 2 (PBS ++ ), fixed at 37 C for 15 min in 4% paraformaldehyde, and permeabilized for 10 min at room temperature with 1% Triton X-100. Cells were blocked with 3% bovine serum albumin, and primary anti-EGFR XP (1:100) antibody was incubated overnight in 1% bovine serum albumin at 4 C. Cells were washed, and secondary anti-rabbit conjugated to Alexa-fluor 488 (1:1000) was incubated 1 h at room temperature. Cells were then stained for DAPI (1 μg/ml) for 5 min and mounted to a microscope slide using Diamond Anti-fade mounting media. After curing, cells were imaged using a 63× 1.4NA oil objective on a Leica SP8 White Light Laser Confocal Microscope. Images were the product of 3-fold line averaging. Three to five images were taken per coverslip, and Pearson's Correlation Coefficient (r) was calculated via the Coloc2 plugin from ImageJ.

Co-immunoprecipitation
A431 cells were plated in a 12-well plate to 80% confluency and allowed to adhere for 24 h. Cells were then starved overnight and treated with serum free-media without (no treatment) or with (PET1) PET1-DL680 at 1, 2, and 5 μM for 1 h before a 5 min treatment with or without EGF (100 ng/ml). Cells were then washed, scraped from the plate using co-IP buffer (50 mM Tris-HCl (pH 7.4), 150 mM NaCl, 5 mM EDTA, and 1% NP40) containing protease and phosphatase inhibitors, and allowed to sit on ice 30 min prior to a 10 min centrifugation at 13,000g. The pellet was discarded and 100 ng of total protein was diluted to 400 μl with anti-EGFR IP (1:100) antibody and rotated at 4 C overnight. Sixty microliters of prewashed Protein A magnetic beads (Cell Signaling) were added and rotated for 2 h. Then, lysate was removed and beads were washed 4 × 10 min at room temperature with co-IP buffer. Protein was eluted at 100 C for 5 min in 2× Laemmli sample buffer containing no dye, and eluate was run on a 4 to 20% SDS-PAGE gel. The gel was imaged for 680 nm fluorescence on an Odyssey CLx Imaging System (LI-COR) before being transferred to 0.2 μm nitrocellulose, blocked with 5% milk, and blotted overnight with anti-EGFR XP antibody (1:1000). The membrane was washed and blotted with IRDye 800CW anti-rabbit secondary and imaged for 680 and 800 nm fluorescence using the Odyssey system as above.
Pulsed interleaved excitation fluorescence cross-correlation spectroscopy PIE-FCCS was used to study the effect of PET1 on the lateral oligomerization of EGFR. Expression vectors from previous studies were used to label EGFR at the C-terminus with EGFP and mCherry (45). These vectors were expressed in COS7 cells purchased from Sigma Aldrich. COS7 cells were cultured in DMEM (Calsson Lab) supplemented with 10% FBS (Sigma Aldrich) and maintained in a humidified incubator with 5% CO 2 at 37 C. To prepare for PIE-FCCS experiments, the cells were split, seeded on to a 35-mm MatTek plate (MatTek Corporation), and incubated until the confluency reached 70%. The plasmid constructs were transiently cotransfected to COS7 cells using Lipofectamine2000 (Invitrogen) approximately 24 h before the data acquisition. Data were recorded on live cells before and after peptide or ligand addition as described previously (45,46). The two-color PIE-FCCS experiment and the auto/cross-correlation analysis allow us to evaluate the expression density, diffusion (reported as an effective diffusion coefficient, D eff ), and the oligomerization state (reported as fraction correlated, f c ). The density of the receptors expressed was in the range of 100 to 2000 receptors/μm 2 . Lower D eff and higher f c values indicate the formation of larger oligomers. To test for the effect of PET1, 2 μM or 2 μg/ml of PET1 were added to the well with 2 ml imaging media and incubated for 10 min before data acquisition. Data were acquired up to 60 min after ligand addition.

Cell migration
A375 cells were plated to 50% confluency on a 10-cm dish and allowed to adhere for 24 h before overnight starvation in serum-free DMEM. Cells were then trypsinized with 0.05% trypsin for the minimal amount of time required to remove from the plate, washed, and brought to a density of 2 × 10 5 cells/ml in serum-free DMEM. Hundred microliter of cells were plated on the top of a 6.5 mm transwell polycarbonate membrane insert with an 8 μm pore size (Corning 3422), while 600 μl of DMEM containing 10% FBS alone (no treatment) or in combination with EGF (100 ng/ml), PET1 (2 μM), or both EGF and PET1 together was in the bottom of the insert. Cells were allowed to migrate 24 h at 37 C and 5% CO 2 before chambers were washed, cells remaining on top of the membrane were scraped off, and cells on the bottom of the membrane were fixed with methanol and stained using H&E. The membrane was then cut from the chamber and mounted to a microscope slide for imaging using a 10× objective on a Biotek Cytation V microplate reader with Gen5 software.

RTK array
A431 cells were blotted using an R&D RTK array kit (ARY001B). Cells were plated at 80% confluency on a 6-well plate, allowed to adhere 24 h, and then starved overnight using serum-free DMEM either alone or containing PET1 (2 μM). Cells were then incubated for 5 min with PET1 alone, EGF alone (100 ng/ml), or PET1 and EGF in combination. Cells were then washed and scraped from the plate using 1X kit lysis buffer and agitated at 4C for 30 min. Lysates were centrifuged 10 min at 13,000g, and the supernatant was quantified using a DC assay kit (Bio-Rad). Two hundred micrograms of total protein was diluted into a total of 1.5 ml Array Buffer and rocked over the preblocked array membrane containing various total RTK antibodies overnight at 4 C. After washing, anti-phospho-tyrosine antibody conjugated to HRP was blotted for 2 h at room temperature, washed, and the HRP was developed using Chemi Reagent Mix (kit). The blots were imaged using the Odyssey system with 10 s to 3 min exposure times.

Modeling of the TM peptides
The NMR structure of EGFR TM dimer (PDB ID: 5LV6) (5) was obtained from www.rcsb.org. The TM region and the membrane proximal N-terminal residues of EGFR from P 641 -I 673 was extracted from the NMR structure, and the remaining modified TM C-terminal residues from C 674 WN 676 were modeled as an extended conformation of amino acids (Φ, ψ = ±120 ) in PyMOL (The PyMOL Molecular Graphics System, Version 2.4. Schrödinger, LLC). The PET1 peptide [T 638 NG PKIPSIATGMVGAELLLGVVAEGIGLFMEEEE 672 ] w a s modeled as TM helix from I 646 -M 668 based on EGFR TM NMR structure as the template using Modeller (47), and the remaining N-terminal residues (T 638 -S 645 ) and C-terminal residues (E 669 -E 672 ) were modeled as an extended conformation (Φ, ψ = ±120 ) in PyMOL. This PET1 peptide has seven mutations as shown in Figure 1A. The SP [F 638 KLAAVNGGIGSTGGIEMVIPGM-PELTLALVEEEE 672 ] was also modeled as TM helix from G 646 -V 668 in a similar manner with the terminal residues as extended conformation.
Inhibition of EGFR via a transmembrane peptide CG molecular dynamics simulation To check the influence of the PET1/SP on the dimerization of EGFR TMs, we built an initial system configuration which contains the two EGFR TM monomers and the PET1/SP peptide, where all three monomer TM units placed perpendicular to the membrane and 5 nm apart from each other. As a reference for the above systems, we also ran simulations with the same system set-up but omitting the PET1/SP peptide, that is, the two EGFR TM peptides by themselves as the control.
The atomistic modeled systems with EGFR TMs: PET1 (2:1); EGFR: SP (2:1), and EGFR TMs alone, each placed as above, were converted to CG representation using the marti-nize2.py workflow module of the MARTINI3 forcefield (48) (version 3.0.4.28) considering the secondary structure DSSP assignment. CG simulations were performed using Gromacs version 2016.5 (49). The setting up of the POPC bilayer was done using the insane.py script (50) (typically 684 lipid and 19,300 CG water molecules for 2:1 peptide systems; 324 lipids, and 4350 CG water molecules for control EGFR TM dimer systems) around the peptides. The pH of the system was 4.5, setting all the Glutamate residues in the peptides to be protonated. The systems were equilibrated for 500 ps. The electrostatic interactions were shifted to zero between 0 and 12 Å, and the Lennard-Jones interactions were shifted to zero between 9 and 12 Å. The V-rescale thermostat was used with a reference temperature of 320 K in combination with a Berendsen barostat at 1 bar reference pressure, with a coupling constant of 1.0 ps and a compressibility of 3.0 × 10 −4 bar −1 . The integration time step was 20 fs, and all the simulations were run in quadruplicate for 4 μs.

Data analysis
Interhelix distances between the COMs of the TM regions were calculated, and the clustering were performed with a cut off 6 Å using the module of the Gromacs by combining the entire trajectories of all the four simulations. PREDDIMER webserver was used for analyzing the TM dimer cluster centers based on the Fscor, helix crossing angle, and helix rotation angle (51). For the latter, the residue L657, which is part of the regular TM helix, was chosen, and angles were read of PRE-DDIMER energy maps. The contact maps for the TM regions between the helices were calculated with a cut off 4 Å for all the backbone and side-chain atoms. Sequence alignments were done using ClustalX (52). Data were plotted in GraphPad Prism (version 6 for Windows, GraphPad Software, www. graphpad.com).
AlphaFold-Multimer (38) was used to predict models for EFGR TMs alone and in complex with the PET1 or SP peptide (2:1 stoichiometries in all cases).

Data availability
Data will be provided upon reasonable request.
Supporting information-This article contains supporting information.